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Abstract 

We consider dynamics of fluctuations in transonically accelerating Bose-Einstein condensates 
and/or luminous fluids (coherent light propagating in a Kerr nonlinear medium) using the hydro- 
dynamic approach. It is known that neglecting the quantum potential (QP) leads to a singular 
behavior of quantum and classical fluctuations in the vicinity of the Mach (sonic) horizon, which 
in turn gives rise to the Hawking radiation. The neglect of QP is well founded at not too small 
distances |x| ^ Ih from the horizon, where Ih is the healing length. Taking the QP into account 
we show that a second characteristic length Ir > Ih exists, such that the linear fluctuation modes 
become regularized for |x| <C /r- At |x| ^ the modes keep their singular behavior, which 
however is influenced by the QP. As a result we find a deviation of the high frequency tail of 
the spectrum of Hawking radiation from Planck's black body radiation distribution, which can 
be described by an effective Hawking temperature decreasing with increasing frequency. Similar 
results hold for the wave propagation in Kerr nonlinear media where the lengths Ih and Ir exist 
due to the nonlinearity. 

PACS numbers: 03.75.Kk, 42.65.Wi, 04.70.Dy 
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I. INTRODUCTION 



Hawking radiation is one of the most impressive phenomena at the intersection of gen- 
eral relativity and quantum field theory. The account of quantum nature of the physical 
vacuum in curved space led to the prediction that a black hole - an object defined clas- 
sically as an object that even light cannot escape - in fact emits radiation,-'^ which is 
characterized by the Planck's distribution of the black body at a certain temperature. 
The derivation involves singular behavior for the eigenmodes of the proper field equation 
in the vicinity of the horizon, which can be normalized only by using assumptions on 
Planck scale physics that cannot be justified completely. For a general discussion the 
reader may consult Ref. Is]. 

Soon after that development it was suggested^ to consider analogous phenomena in 
condensed matter physics where the "high-energy" (short-wavelength) physics is known. 
The suggestion was based on the observation that the derivation of the Hawking radiation 
uses only the linear wave equation in curved space-time, and does not make the actual 
use of the Einstein equations for gravity. The same conditions for wave propagation arise 
at the consideration of sound propagation in the fluid when the background flow is non- 
trivial.- In particular, the background may be that of a stationary accelerating transonic 
flow where a surface appears, on which the fluid velocity equals that of sound, i.e where 
the Mach number is one, M = 1. This surface has many features similar to that of the 
black hole horizon and is often called sonic or Mach horizon. 

Currently theoretical studies of such artificial black holes became a very active field 
involving a variety of physical systems.-"— Experimentally a white-hole horizon was ob- 
served in optical fibers,— where the probe light was back-reflected from a moving soliton, 
and a black-hole horizon was observed in a Bose-Einstein condensate (BEC) system.— 
Continuing the study of optical fibers a radiation has been quite recently observed,— 
which is a very promising contender to being analogous to Hawking radiation in a table- 
top experiment. A "horizon physics" is currently intensively studied also in the surface 
water waves.—"— 

Creating a stationary laminar transonic flow of an ordinary fluid, as required by the 
original proposal by Unruh,"^ may be a formidable task that is why many authors turn 
their attention to the BEC (see, e.g. Refs. IsMlOl). or as recently proposed, to a "luminous 
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FIG. 1: A scheme of Laval nozzle as a vessel with non-monotonously varying cross-section. It 
is capable of accelerating a flow from subsonic to supersonic velocity. The sound velocity is 
reached in the throat (the narrowest cross-section) as shown by the dashed line where a Mach 
(sonic) horizon is formed. 



yi2oi) schematically shown in Fig. [1] accelerates a flow from sub- to supersonic velocities. 
It uses the fact that the flow accelerates if the cross-section of the vessel decreases in the 
subsonic region, whereas in the supersonic region the cross-section must increase in order 
to accelerate the flow. An optical analog of the Laval nozzle^^ is a properly shaped wave 
guide fllled with a Kerr nonlinear medium, in which a coherent laser light propagates. 

BEC is usually described by the Gross-Pitaevskii equation, which can be converted 
into two equivalent hydrodynamic equations by means of the Madelung transformation.— 
A similar procedure exists also for the coherent light.— Analysis of the behavior of fluctua- 
tions near the Mach horizon based on these equations is usually carried out neglecting the 
role of the quantum potential (QP).-'^^ This approximation can be well justified nearly 
everywhere except for the narrow vicinity of the sonic horizon due to the singular behav- 
ior of the fluctuations in this region. Neglecting QP the fluctuations are described by 




fluid"i^ 



in an optical analog of the Laval nozzle. Generally the Laval nozzle (see, e.g., Refs. 
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an equivalent Klein-Gordon equation for the phase fluctuations in a curved space with 
Schwarzschild metric, as flrst proposed in Ref. \^. The reader should also note that the 
fluctuation modes are then dispersionless, i.e. their frequency depends linearly on their 
wave-number, if QP is neglected. 

At small distances from the horizon the contribution of the QP (kinetic energy oper- 
ator for atoms in BEG or dispersion in luminous fluid) diverges and cannot be anymore 
neglected. This reminds the behavior of vacuum fluctuations near the black hole horizon 
on the Planck scale. Since in our case the physics on small scales is known the role of the 
QP in the near vicinity of the Mach horizon can and should be addressed. 

Accounting for the QP introduces dispersion effects at higher frequencies corresponding 
to small length scales. The question how far such effects influence Hawking radiation has 
been studied during the last two decades. It seems that the flrst such study was done 
by Jacobson2^i2i. Whereas Ref. |23| had stated that a derivation of black hole radiation. 



taking the ultrashort distance physics into account, was missing, such a derivation was 



presented in Ref. ]2S- It has been concluded that boundary conditions can be formulated, 
which lead to the Hawking effect. In order to deduce quantitative results, e.g. for the 
flux and the spectrum, the mode equations based on speciflc dispersion laws were solved 
numerically^>2^"— and analytically.—"— In all of these papers clear evidence has been given 
that Hawking radiation is basically a low frequency phenomenon and is robust, although 
counterexamples were presented^ predicting strong deviation from the Hawking's result. 
Due to dispersion an upper threshold comax exists in the superluminal case such that 
the Hawking radiation disappears for cu > cumax-—^ For more details the reader is also 
referred to the review in Ref. 38. 



In the current paper we will focus on BEG or equivalently on the light in a Kerr non- 
linear medium. Assuming a stationary background solution with inhomogeneous density, 
p(r), and velocity, v(r), flelds we will derive the linearized equation for the fluctuations 
of both flelds, 6p{r) and 5v(r), accounting for the QP. It will have a form more compli- 
cated than the Klein-Gordon equation of Ref. \4- Reducing the problem to a spatially 
one dimensional one, two Hermitian flelds xi^^t) ^^id C,{x,t) are introduced, which obey 
two partial differential equations. We show that assuming linear dependence of the back- 
ground velocity and density on the coordinate near the horizon x = 0, the equations can 
be decoupled in the momentum space, and reduced to an equation for Xk representing the 
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density fluctuations Sp{r). Our mathematical analysis is close to that of Ref. |3ll, which 
makes use of the Laplace back transform. However, there are several formal and physical 
points of difference from Ref. 131 



which will be discussed below. In addition the resulting 
dispersion relation is dictated by the GP (NLS) equation ([T]) (see next section) with the 
coordinate dependent background flow rather than chosen by the mathematical require- 
ment of nonlinearity^^"— >^""— or assuming a rather general nonlinear dispersion.— i^i^^ 



II. EQUATIONS FOR BEC AND LUMINOUS FLUIDS 

We start by presenting a brief overview of the basic equations controlling the dynam- 
ics of BEC or light in Kerr medium. This will allow us to make necessary definitions 
and introduce basic notions needed for the further analysis. Dynamics of Bose-Einstein 
condensate is quite successfully described by the Gross - Pitaevskii (GP) equation 

ihdM^, t) = V2$(r, t) + f/e.t(r)$(r, t) + g^r, t)\^^v, t) (1) 

where $(r, t) is the condensate wave function of identical Bose particles with the mass m, 
Uexti^^) is an external potential, and g is the parameter of the particle interaction. The 
equation is usually derived (see, e.g. Ref. 13) from the equation of motion for the boson 
field operator \E'(r,t), which is represented in the form 

$(r,t) = <l>(r,t)+V^(r,t) (2) 

so that the field operator il){v,t) describes fluctuations near the condensate state $(r). 
For small fluctuations the corresponding linearized equation reads 



ihdt^iv, t) = ("I^V^ + f/e.t(r) - /i^ i^(r, t) + 2g\^{v)\^^{v, t) + g^\v)^\v, t) (3) 

Usually solution $(r) of the stationary GP equation is used, which is obtained from Eq. 
([1]) by replacing its left hand side by /i$(r), where /i is the chemical potential. 

A similar mathematical description is applied to the coherent light propagating in Kerr 
nonlinear medium (see, e.g., Ref. |40[ ). The corresponding equation 



id,A 



(4) 



is usually called Nonlinear Schrodinger (NLS) equation. The latter is derived from the 
classical Maxwell equations for electromagnetic waves in a Kerr nonlinear medium. The 
paraxial approximation is assumed, according to which the electric component of the wave 
is represented in the form 

where the amplitude A weakly depends on the coordinate z along the light propagation 
axis. As for its frequency dependence it is narrowly peaked around the principal frequency 
Uq. Here f3o = k{uo) where k'^{uj) = u^fioe^u). The magnetic susceptibility /io is assumed 
to be a constant, whereas the frequency dependence of the dielectric function e{u) plays 
an important role, c is the light velocity. Usually the light polarization does not change, 
which allows one to treat this equation as a scalar one. 

Although NLS Eq.(j4]) is classical we can formally multiply it by h and get the corre- 
spondence to the quantum GP equation where the propagation distance z plays the part 
of time, hPo is the mass of a fictitious particle, which is measured in units of momentum. 
It is interesting to note that a quantity corresponding to energy or chemical potential 
should be also measured in units of momentum. It implies also that the "velocity" is 
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dimensionless. Spacial variation Uext{x,y) = —h^^^An'^{x,y) of the refraction index 
= fio£{uj)c'^ plays now the part analogous to that of the external potential in the GP 
equation and the field dependence of the dielectric function results in the nonlinear term 
in Eq.©. 

The time coordinate appears now as one of the "space" coordinates in the effective 
Laplacian 

V' = dl + dl±d'^ (5) 

where 



1 ^ dk{u}) 



is the light group velocity in the medium and D = ^ d_k(uj) 



is the 

UJ=UJo 

dispersion coefficient. The minus sign in Eq. ([5]) corresponds to the normal dispersion 



2tt dui^ 

LJ=UJO 



{D > 0) whereas the plus sign corresponds to the anomalous dispersion {D < 0). The 
hyperbolic "Laplacian" for the normal dispersion can often occur for the coherent light 
propagating in a Kerr medium but does not appear in the GP equation for BEG gases. 



We can always choose the frequency window within the anomalous dispersion regime such 
that the Laplacian keeps its standard shape for the coherent light as well. When a laser 
emits a pulse, z is actually the position of the pulse reached during the time zjvg. t 
is used to describe the shape of the pulse (in the z direction) in the coordinate system 
moving together with it. 

The amplitude A can be always rescaled in such a way as to make its dimension 
coincide with that of the wave function $. It does not influence the linear part of Eq.(|l]), 
whereas the factor appearing in the nonlinear term can be absorbed into the parameter 
A, which together with h then maps onto g of Eq.([T]), making the mapping between the 
two equations nearly complete. 

Now we can also consider stationary solution of the NLS equation and small fluctuations 
around it, 

A{r,z) = ^{^ +i;{r,z). (6) 

where r = {x,y,T). tp(r,z) satisfies an equation, which maps onto Eq. according 
to the above prescriptions with the only difference that ip in Eq.([6]) is a scalar function 
rather than a field operator of Eq.([2]). 



III. EQUATION OF MOTION FOR FLUCTUATIONS 

The hydrodynamic description of the Heisenberg equation of motion for the field op- 
erators ^l/(r, t) and \&^(r, t) follows by introducing the density operator 

p(r,t) = m^\r,t)^{r,t) 

and current density operator 

J(r,t) = i[$t(r,t)V§(r,t) - (V§t(r, t))$(r, t)]. 

Both are connected by the continuity equation 

dtp{r,t) + V}{r,t) = 0. 

Substituting Eq.([2]) into Heisenberg equations of motion for '^{r,t) and \E'"''(r, t), us- 
ing $(r,t) = f{r,t) exp[—iip{r,t)] and neglecting fluctuations one obtains hydrodynamic 
equations 

|p + V ■ [pv] = 0, (7) 



|v + ivv^ = -i-V %u + + gf] (8) 

for the density p{r,t) = mp{v^t) and velocity v(r,t) = — ;^Vv?(r, t); both / and are 
real functions. These two equations are equivalent to the GP equation ([1]). Here 

is QP. This hydrodynamic representation of the basic equations allows us to think about 
the BEG fluid or luminous fluid, respectively. 

We now consider small fluctuations of these two fields with respect to a stationary 
solution $o(i") = /o(r) exp[— 'i(y9o(r)]. It follows from Eq.(I7]) that 

V ■ [po(r)vo(r)] = 0, (10) 

where po{r) = mfQ^r) and vo(r) = —^Vipo^r). To describe the fluctuations we introduce 
the dimensionless quantities 

X = ^[e-^°^^ + e*^o^], ^^^^ 

which are field operators in the case of BEG or classical scalar fields in the case of coherent 
light. Therefore we skipped the hats. These quantities are connected with the density 
and velocity fluctuations by the relations 

Mr,t) = --Ve(r,t), 

m (12) 

Sp{r,t) = po{r)x{r,t). 
Then it is straightforward to show that Eq.(|2]) and its Hermitian conjugate lead to 

5x--^v(/oVO = o, 

(13) 

where D = dt + vo ■ 'V . These equations also follow from Eqs.© and ([8]) in the linear 
order in 6p and 5v. While deriving equations ( IT3|l . we used Eq. lfTOj) . 



It is instructive to indicate 
discussion in Appendix of Ref. 



;hat Eqs. ( fT3l) are Euler - Lagrange equations (see also 



28l ) corresponding to the Lagrangian density 



^ = Ifoixi - ex) + ^/o^o ■ (xVe - eVx) - - ^/o (VO' - ^/o (VX)' (14) 



This Lagrangian density allows us also to define canonical momenta 



Px 
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5? 



/oX (15) 



for each of the fields x ^"^^ ^- Together they form two canonical pairs, each of which 
can be used to derive Hamiltonian. The Hamiltonian density for these two cases can be 
derived as 



n 

1 



7oxe 



c 



foX^ 

h 



^/o^vo ■ (xve - evx) + ^^?/oV + ^/o (vo^ + ^/o (vx)^. 



(16) 



Although the derivation (flGll depends on the choice of the canonical pair, the final Hamil- 
tonian does not depend on it. Eqs. f|T3|) are now just the pair of Hamilton equations 
of motion for the Hamiltonian density (fT6|) . The above manipulations with adding or 
subtracting a full derivative indicate that the variations of the Lagrangian with respect to 
X and ^ are not uniquely defined. This feature of the Lagrangian (HM does not however 
effect the Euler- Lagrange equations f lT^ but is important when defining the canonical 
momentum and Hamiltonian. It means that we have to choose which of the functions x 
or ^ plays the part of the canonical coordinate, then the other one will enter the definition 
of the conjugate canonical momentum. It is emphasized that we can use only one of these 
pairs at a time, but not both of them simultaneously. The physical result is independent 
of the choice. 

If ^ is chosen as the canonical coordinate then the quantization condition reads 



[p^(r,t),e(r',t)] = [f^xir,t),ar',t)] = M(r - r') 



(17) 



where [■■■,■■■] denotes the commutator of two operators. The other choice of the canon- 
ical pair would obviously yield the same quantization condition. One can readily see that 
the quantization condition ( |T71) in the BEG case is consistent with the definition ( fTTj) . 
where ip"^ and ip are Bose field operators. In case of the coherent light when these are 
scalar functions, Eq. ( !T7|) allows us to quantize them thus removing the last remaining 
difference between BEG and luminous liquid. 



We may extend the definition of tlie fields x ^i^nd ^ and consider them complex. Then 
the Lagrangian density f|T4l) can be represented in the form 



where the two component field 

d= I ^^'^ I (19) 




has been introduced. 

The corresponding equation of motion reads 

^/o^.^ + ^/o vo ■ cjyVd - j^gfo\l + a^)^ + ^V(/oV^) = 0. (20) 

It is obviously equivalent to two Eqs. (fT3l) . Now we multiply Eq.( |20l) by 2??^ from the left, 
write also a Hermitian conjugate equation and take their difference. The result is 

5t(/o^^M) + V(/oVo^V,^) - ^V[fl[d^Vd - = 0. 

We may conclude that the quantity q = f'^d'^ayd plays the role of the density of this two 
component field, whereas 

j = /o vo^V,^ - ^ A/2[^tv^ - (V^9t)^] (21) 

is the corresponding current. The integral 

d\fl§^ayd = -i j d\fl{x*i - Cx) (22) 

is a conserved quantity, which can be used in order to normalize solutions (eigenfunctions) 
of Eq.f llSp or to define a scalar product of two such solutions. Similarly to the well-known 
Klein-Gordon norm, Eq.( l22|) is not positively defined. Positive and negative frequencies 
are considered separately in order to resolve this problem in the case of the Klein - Gordon 
equation (see, e.g. Ref. [si). A similar procedure will be applied in this well. This 

property will be important for the further discussion. 
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IV. EIGENFUNCTIONS 



A. Bogoliubov spectrum 

As an example we first may consider a homogeneous background solution $(r), in which 
only the phase may linearly depend on the coordinate, i.e. both the density po = ^fo 
and the velocity mvo = —KVifo are constants. Then by solving the second equation in 
Eq.( lT3|) for x and substituting the result in the first one this system reduces to a linear 
partial differential equation with constant coefficients: 



D'^ - —,PoV'i + ^ {^y') \ = 0. (23) 



It straightforwardly leads to the celebrated Bogoliubov spectrum 



[u; - vo ■ k]' = s^k^ 



1 + ^ 



n\k) (24) 



in the coordinate system moving with the constant velocity vq. The quantity s ^ = ^ is 
obviously the sound velocity in the long wave limit and = fv^ / {2mgfQ) = h^/ {2m?'s ^) is 
the healing length in BEC and the nonlinearity length = 1/ (2/3oA|Ap) in optics. These 
lengths are determined by comparing the QP Uqu with the nonlinear term in the equation 
of motion. Eq. fl2^ demonstrates that the propagation at high frequencies corresponds 
to the so called superluminal case. For Ih = 0, i.e. neglecting QP, Q{k) reduces to a linear 
dispersion law. It is interesting to indicate that the Bogoliubov spectrum of excitations 
appears not only in the quantum BEC system, for which it has been originally proposed, 
but also in the case of classical electromagnetic wave propagating in the Kerr medium. 
Eq. flMl) has four solutions, which for vq = read 



k' = -^ l±\l + ^] (25) 




Two of the four solutions (— sign in Eq.( l25|) ) are real and correspond to two plane waves 
moving in the mutually opposite directions. The other two solutions (+ sign in Eq. fl25|) ) 
are imaginary and correspond to excitations growing or decaying exponentially on the 
scale Ih- They are usually of a minor physical importance. The situation may change for 
the supersonic regime when four real solutions are possible, provided u < Ur, 



27^ 
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B. Hydrodynamic approximation 



The hydrodynamic approximation, apphed in many cases, assumes that the QP Uqu in 
Eq.([H]) is neglected. It is justified when the QP is small as compared to the interaction 
energy, Ugu ^ gfo, i-e. on the large scales |x| ^ Ih- Now again it is possible to eliminate 
the variable x iii Eq. (fT3|) . which is then reduced to the equation 

^{[d^ + dtvo ■ V]e + V(-vo9iO + V(-vo(vo ■ V))e - V(s2V)e} = 0. (26) 
5' Jo 

for the phase fluctuation fleld C,- 

This equation has the form of a Klein-Gordon equation in a curved space, which has 
been discussed many times starting from the seminal paper by Unruh^ (see also Refs. 

In particular, we may apply it to the situation, which takes place in the 
Laval nozzle (Fig. [1]), creating a transonic flow of a luminous fluid.— In this case we 
may reduce the problem to 1+1 dimensions keeping the time coordinate and one spatial 
coordinate x along the streamline crossing the Mach horizon at x = 0. In its narrow 
vicinity the velocity is a linear function of the coordinate, Vo{x) ~ s{l + ax). Since the 
density po{x)vo{x) = ^s'^{x)vo{x) of the flow reaches its maximum at the throat of the 
nozzle we have s{x) ~ s(l — Therefore the functions v{x) and s{x) cannot be chosen 

independently unless we allow for the x dependence of the nonlinearity coefficient.— Then 
Eq.( l26|) becomes 

{dl + 2(1 + ax)drd^ + adr + 3ad^{xd^)} ^ = (27) 

to within the linear terms in x, where the variable t = ts is used. This equation has two 
solutions, which at — )■ behave as 

6 oc e^°'"^-^'^*, ^2 oc e^'^^-^'^* (28) 

where 70 = and k{ijj) = f 2w^M^' Correspondingly we have also two eigenf unctions 
for the density fluctuations 

Xi oc e(^«-i)i'^"-*^*, X2 oc e''^"-'"* (29) 

These two solutionai^i^^ are the former propagating plane waves in the Bogoliubov 
spectrum strongly distorted by the accelerating background flow. The solution ^1 corre- 
sponds to the excitation "attempting" to propagate against the flow in the vicinity of the 
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Mach horizon, which cuts it into two parts moving in the opposite directions from the 
horizon. The solution ,^2 is a fluctuation propagating with the flow and therefore have 
a shape quite similar to that of a plane wave. In the high frequency limit these modes 
propagate with the double sound velocity. 

Two remarks are in order: The exponent 70 here differs from the corresponding ex- 
ponent in Ref. [ij by a factor of two due to the misprint in that paper. This exponent, 
as well as the Hawking temperature to be obtained below, differs also from the exponent 
in Ref. 
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3J by a numerical factor due to the fact that here we take into account the 



coordinate dependence both of the velocity of the flow and of the local sound velocity. 

An important feature of these two solutions is that they behave completely differently 
at I a; I — )■ 0. The solution ^2 is a smooth function in this region meaning that the corrections 
due to the neglected QP, proportional to the second derivative of ^2 are small and the 
above approximation is quite sufficient. Completely different situation takes place when 
dealing with the solution ^1, which is singular in the limit \x\ — )■ 0. The corrections to ^1 
due to the QP diverge at |x| — )■ and a more refined procedure is necessary in order to 
regularize the behavior of ^1 in this limit. 



C. Regular izat ion 

As indicated above the eigenfunctions ^1 and xi, obtained in the hydrodynamic ap- 
proximation, (1281) and (129|) . have a singularity (branching point) on the horizon at x = 0. 
This type of behavior is typical of the eigenfunctions of the Klein-Gordon equation in the 



vicinity of the horizon of a real black hole (see, e.g., Ref. |43| and references therein) and 
leads to the celebrated Hawking radiation. Regularization in this case requires knowing 
the physics at the Planck length scale (see, e.g. Ref. Q). In the systems discussed here the 
part of Planck length is played by Ih, which is the healing length for BEC or nonlinearity 
length for the luminous liquid. In order to regularize Xi? near the horizon we have to 
take into account the QP in Eq. (fT3|l . 

We apply now the same assumptions concerning the coordinate dependence of the 
velocity and density near the Laval nozzle throat as in the previous subsection. Using the 
Laplace representation 

X{x,r) = e-'^^ [ dkxke'"'' 
'c 
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where v = u/s^ and the contour C is chosen such that the integral converges,— then Eqs. 
( 1T3|) become 



h 



h 



-iu + isk)xk - asdk{kxk) + iak—^k H k = 



m 



m 



(30) 



[iak + k^] Xk - ms^iXk - iadkXk) + h[{-iuj + isk)ik - asdk{k^k)] = 0. 



4m 

where the terms 0{a'^/k^) have been omitted. This approximation is equivalent to consid- 
ering a region a\x\ <^ 1 near the horizon. The representation fl3Ql) is especially convenient. 
It allows us to solve explicitly the first equation, 

m {iu - isk)xk + asdk{kxk] 



h iak + k"^ 

and substitute the result into the second equation. We then get equation 



dk In Xk 



-i{iak + ef-e + {v-kf 



(31) 



(32) 



iak{2i' — 3k — ia) 

for the function Xk- Two comments are in order with respect to the approach taken in 



Ref. 



31 



where the Laplace representation has also been used. Although Eq. (1321) has the 



structure similar to that of Eq. (23) in Ref. 1311, it is certainly different. Most important 



31 



uses 



is that the r.h.s. of Eq. |32]has two poles instead of a single pole at A; = (Ref. 
k = —is). This difference stems from neglecting one of the terms in equation (12) in Ref. 
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However, as we will see below, it is crucial for our analysis not to neglect this term. 
An attention is also drawn to the fact that the r.h.s. of Eq. [32] involves a cubic term with 
purely imaginary coefficient. 

Eq. fl32l) can be readily solved to within a term independent of k, 



2 i 

In Xk = 7i + 72 ln(A; - - 3") + H^, ^) 



(33) 



where 

1 iu 
^^ = 4-2^' 
and the Ih dependent part is given by 



72 



1 -1 4x 2 s 14,2 2 

t — u Itu^ H Itu^ 

4 6a 81a ^ 81 ^ 



Afk, I/) = S + —ak^ - — z/P - —u^k + —uak 

^ ' ^ a 1 18 36 18 27 27 



(34) 



The real space solution of Eq. f[T^ is finally obtained by the back Laplace transform 



dkk 
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c 



k u a 

3 3 
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exp {A(A;, v) + ikx} . 



(35) 
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This integral is ca. 



functions in Ref. 



41 



culated in the way similar to that used for the calculation of Airy 
First, we have to find the sectors of convergency of this integral in the 
Riemann plane for the complex variable k, which are controlled by the term proportional 
to k^ in Eq. flMl) . Hence the condition ImA;'^ < must be required. It holds in the three 
sectors: 2tt/3 > arg/c > 7r/3, Att/S > arg/c > tt and 2tt > axgk > Stt/S, shown in Figure 
[21 The integrand's function has also two branching points at A; = and k = 2z//3 + 2a/3, 
which result in two cuts shown in Figure [2l Due to these two cuts there are five possible 
integration contours Ci to C5 and any other contour may be represented as a combination 
of these five. They provide us with five different results for the integral (1351) . however the 
sum of these five integrals is zero due to the Cauchy theorem. As already mentioned 
above a similar approach to the problem was applied in Ref. |31| with nearly the same 
configuration of contours as in Fig. [2j However, the corresponding function in Ref. l3l| has 
only one branching point. But to have four independent solutions we need two branching 
points rather than one. They are also important for Eq. ( !37l) to be derived below. 

The corresponding four integrals may be calculated by the steepest decent method. 
Then each of these integrals will be a linear combination of the contributions due to the 
four saddle points, determined by the equation 

dklnxk = -ix. (36) 

Alternatively we may also choose the integrals due to the individual saddle points as 
linearly independent solutions. Substituting In Xk from Eq. fl32l) into Eq. fl36|) we get 
after some manipulations the x dependent dispersion relation 

+ A;') V + k^s^ix) = [u - kvix)f = fi'(A;, u) (37) 

where we remind that v{x) ~ s(l + ax), s{x) ~ s(l — \ax^ and uj = su. Here some terms 
of the order 0{{a/kY) and 0{{axY) are kept only for the sake of a compact presentation. 



■ul. 

and Eq. (23) in Ref. 



31 



28. 



However only the terms up to the first order are meaning. 

Equation fl37l) is rather similar to Eq. (6) in Ref. 
Nevertheless we have to indicate important differences. Those equations relate to the 
regions a\x\ ^ 1 far from the horizon where it is assumed that the flow velocity reaches 
its supersonic or subsonic limits, and the corresponding velocities, f=p, do not depend 
on the coordinate x. Our equation (!37j) is deduced for the region close to the horizon 
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FIG. 2: Five possible contours in the complex k plane for calculation of the integral ()35p . The 
integral diverges in the dashed sectors. The contours lie in the sectors of convergency and also 
go around the two cuts, which start from the two branching points at k = and k = + |a. 

a\x\ < and both the flow and sound velocities depend on the coordinate. It also 
contains an imaginary correction in the quartic term. 

An interesting insight may be obtained from this equation by a graphical analysis, 



which is similar in some aspects to the one carried out in Refs. |30|J31| and particularly in 
Ref. |28|. Let us neglect the imaginary correction to the quartic term, which can be taken 
into account as a perturbation, and plot both branches of the function Q{k,uj) crossed 
by the three straight lines u — kv{x) for three different values of Fig. [31 We also 
do not take into account the x dependence of the sound velocity, since it would have 
overloaded the graphical presentation without changing any conclusions. All the three 
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FIG. 3: (Color online) Graphical analysis of equation (I37p . Two branches of the function Q{k,uj) 
are plotted in red in figures (a) and (b). They are crossed by straight dashed lines uj — kv{x) for 
three values of x < ((a) - subsonic regime) and three values of x > ((b) - supersonic regime). 
The full circles show the most important solutions of the equation. The long dashed lines show 
the function Q{k,uj) in the absence of the dispersion, i.e. = 0. 

lines are plotted for a fixed frequency but decreasing distance |x| from the horizon on 
both the subsonic side (Fig. |3^) and the supersonic side (Fig. |3]d). The line (c) in both 
cases corresponds to |x| =0 when v{x) — > s. 

In the subsonic region (Fig. ISti) where x < we always have two real solutions r and 
s, with and without dispersion. It is interesting to indicate that in the absence of the 
dispersion, i.e in the absence of the quartic term in Eq. (!37|) when ^ and when the 
two branches of Q{k,u) are just two straight lines with the slope one in dimensionless 
units (long dashed lines in Fig. [3]), the s solution moves to infinity {kg — oo) as x — > 0. 
The dispersion changes this behavior so that kg tends to a finite value in this limit (going 
from line (a) to (b) and then to (c)) at the point (s') where the (c) straight line crosses 
the Q{k, (jj) curve. 

In the supersonic case, the straight line (a) crosses the curve Q{k,u) in four points 
producing three right movers and one left mover (not all are visible in Fig. [3]). It is 
important to emphasize that two of these solutions, i.e. both right movers (r) and (s) 
exist even in the absence of the dispersion. If we move closer to the horizon, |x| — )■ 0, the 
slope of the straight lines tends to that corresponding to the sound velocity (line(c)). In 
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the absence of the quartic term in Eq. (l37Il the s solution moves to infinity as we approach 
the horizon a; — )■ 0. 

Two solutions (both for subsonic and supersonic case) 



Sax' ' 4 — a 

are obtained directly from Eq. ( 1371) at Ih = 0. The saddle point kr leads us straightfor- 
wardly to the right mover eigenfunction X2 of Eq.f l29l) . which will be now denoted by Xr- 
As for the kg point, the integral fl35l) converges in a very broad range of k values. We may 
introduce the new variable z = kx and write the integral (l35l) in the approximate form 

(x, r) ^ a;~(^^+^^+i)e-*^" / dzz^^^^^ exp {iz} . (38) 



Xs 

where neglecting A(fc, v) in the exponent of the integrand is justified as long as x S> Ih- 
(A stricter inequality will be given below). The validity of this equation is limited also 
by the inequality x <^ min{l/z/, As a result we get the singular eigenfunction 

7-1 

in which the parameter 

7 = -71 - 72 = ^ + - 39 

Sa 81 a 27 

differs from the parameter 70 in Eqs. (1281) and fl29|) for ^1 and xi, respectively, since it 
contains small corrections due to the QP. That is why the notation Xs is used instead of 
Xi- This type of singular behavior of the eigenfunction typically occurs near the horizon 
and is a crucial ingredient in formation of the Hawking radiation ( see, e.g., Refs. I42|j43l ) . 

Accounting for the quartic term in Eq. ( 1371) introduces important changes. While 
moving towards the horizon, |a;| — > 0, for a given frequency we rotate the straight lines in 
Fig. |3]so that the line (a) becomes (6) when two real solutions merge (shown as the point 
s' in Fig. [3]b), and then disappear (become complex) (c). It actually means that the 
singular behavior (138 p of the eigenfunction Xs holds for not two small values of |x| and 
becomes regularized closer to the horizon. Here a remark is in order. It was indicated in 



Refs. 



27 



28 



that there is a maximal frequency Umax, above which this pair of solutions 
does not appear at all, which makes the Hawking radiation impossible. This condition 
was obtained by considering the behavior of the eigenfunctions in the asymptotic region 

18 



X — )■ oo far from the horizon where the flow velocity becomes constant. Our equation (|3j 
allows for a smooth interpolation between these two regions. It is clear that if we start at 
uj > Umax in the supersonic asymptotic region those two solutions do not show up, and 
then moving closer to the horizon, leading to a variation of the sound and flow velocities, 
the singular solution s will have no chance to appear in the close vicinity of the horizon 
as well. If in contrast u < Umax, both solutions exist for x oo. Now approaching 
the horizon a critical distance x^,{u) exists such that that these solutions disappear for 
X < x^:{co). 

We may put the above pattern on a more quantitative basis. In order to find the 
behavior of Xs in the nearest vicinity of the horizon we substitute kg ^ x~^ in the dimen- 
sionless exponent A in Eq.f lS^ and obtain that the cubic term in Eq. flM|) (former quartic 
term in Eq.f l32p ) becomes large when 

|x| < = (40) 

In this region the convergency of the integral f l35|) is controlled by the cubic term in A and 
the result remains finite and not singular at |x| — )■ 0, i.e. the function Xs is regularized at 
I a; I < Ir- In order to demonstrate this it is sufficient to calculate the integral 

f°° 1 

i / 2;~'''exp( l^z"^ — zx)dz = 

Jo 18 

3 fi 




3 J \ 3 ' 3'3' 3/3 
whose integrand is the approximate integrand of (1351) at large k. It represents the integral 
(135|) over the upper half of the contour Ci or C5. Here 1-^2(0; b, c; x) is a hypergeometric 
function and r(a;) is the Euler Gamma-function. The other integrals corresponding to 
other contours or their parts can be obtained by changing the relevant phases. This 
integral tends to the finite limit 

1^ !|V^expf-lo.Vfi^l (42) 



3 V 18 y v 2 
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for X — !■ 0, in fact for |x| ^ Ir- The derivatives with respect to x in this hmit can be 
also calculated, resulting in a regular Taylor expansion. One can also readily see that this 
function may have zeros at complex x. 

It is emphasized that the integral fKT]) depends on the regularization length rather 
than on the healing length 1^. It is also interesting that typically lhC( < 1 so that regular- 
ization starts at the scale larger than the healing length, 1^ > Ih- 



V. HAWKING RADIATION FROM THE MACH HORIZON. 

In order to calculate the frequency spectrum of the Hawking radiation emanating from 
the Mach horizon we will apply the approach similar to that proposed by Damour etal.— 
The central point in this approach is the calculation of the norm of a straddled fluctuation, 
which clearly demonstrates how a negative frequency state is cut into negative and positive 
frequency states propagating in the opposite directions from the horizon. In our case we 
have to calculate the norm (122]) for a system with variation in the x direction only: 

< >= j dxQs{x) = j dxfl{x*,is - CsXs) (43) 

for the pair of eigenfunctions {S,siXs)- We have to deal here with the two component 
function (I19p . of which Xs calculated above is only one component. The second component 
can be found now, say, by solving Eq. ( |T3l) . Neglecting again the contribution of the 
QP in Eq. (II3D we get 



Xs = TTlDis (44) 



so that 



6 « X 



7 



(45) 



for min{l/z/, !/«} ^ |x| ^ Ir- 

The two component field density and current density are 

Qs = -^{mes)^s - &^s] + t^o(x)[(o'.C)6 - CsdM} , 
Js = t {voi&^s - (5tC)e.) + (vlix) - s\x)md.^s - (d^O^s)} . 

where the factors h and g are absorbed in the normalization factor of C,s- The reader 
should also note that a second order term of Xs, which is 0(/^) has been neglected in js, 
to be consistent with the neglect of Ugu- One can readily see that the definition of the 
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norm (j^3|) coincides in this approximation with the Klein- Gordon scalar product in the 
corresponding curved space. 

The local coordinate transformation 

,~ Vo(x)dx ,~ , 

at = at -\ — 9^T; "'^ — (46) 

s^{x) — Vq{x) 

proposed in Ref. |j can be also written as 

dx''= [ I = A^.dx" 
dx 

It allows us to rewrite the density and the current in the form 

^ = - %^a'^ga l'^-'^--g^''- (47) 
3, = ~i[vUx) - sHxMd^C.K. - QdiU 

which have different signs in the subsonic v(x) < s(x) and supersonic v(x) > s(x) regions. 

Now we use the density (H7|) in Eq. ( H3|) in order to calculate the norm. The corre- 
sponding integral is separated in two regions outside the "black hole" x < —Ir (left) and 
inside it x > (right) (see Figured]). We may neglect the contribution of the narrow 
region |x| < in the integral where full solution of the second Eq. ( lT3l) . rather than Eq. 
OH]), should be used. As a result the norm takes the form 

— Ir /"OO 

dxQs+ / dxQs=<'&s,'&s>left + <'&s,'&s>right ■ (48) 

The function 'qs in Eq. fHT]) diverges at \x\ — )■ 0, which could have resulted in the divergent 
contribution of the neglected region of integration. However, it is important to emphasize 
that Eq. ( H71) holds only outside the narrow region \x\ > Ir near the Mach horizon. Within 
this region, \x\ < Ir, we have to return to Eq. fH3]) . Although infinities are indicated as 
the integration limits in Eq. ( HH]) . the principal contributions come from the regions 
min{l/a, l/i^} > \x\ > Ir- Since both Xs and C,s are regular at |a;| < Ir we will get only 
a small correction to the norm fl47p from this narrow region, which can be neglected as 
long as Ih <^ min{l/a, l/i^}. 

The integrals on the two sides of the horizon (left and right) approximately obey the 
relation 

< ^s,^s >ieft= -e2^I"^7 < >^^^^^ 
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due to the analytical properties of the function The total negative norm ( HHj) is cut 
into the left moving positive frequency state and the right moving negative frequency 
state. The left state propagates against the flow "outside the black hole". Its relative 
weight is 

where 

3/3 



Im7 = — — h 



2u 1 uH 



r 



3sa 162 s ^ 

contains a correction due to the QP proportional to the third power of the frequency and, 
hence, N{u) deviates from the Planck distribution for the black body radiation. This 
result can be also understood as a black body radiation spectrum with the frequency 
dependent temperature 



1 + TFT5^ 



108 s'^ 



For not very high frequencies u <^ s/l^ this dependence may be neglected and A^(co') 
becomes the standard Planck distribution with the effective Hawking temperature 

?>tisa 



Th{0) 



VI. CONCLUSIONS 

The paper analyzes the role of the QP in the behavior of straddled fluctuations in 
the nearest vicinity of the Mach horizon. These are fluctuations attempting to propagate 
against the transonic flow that rips them into two parts, one of which propagates against 
the subsonic flow outside the "black hole", whereas the other part is taken away by the 
supersonic flow and "falls down" inside the black hole. The fluctuations are characterized 
by the singular behavior fl28l) in the region ^ a; <^ min{s/a;, which on the 

quantum level leads to the Hawking radiation from the horizon. The singularity at the 
horizon |x| — )■ is regularized due to the QP at the scale x < Ir (I^Oj) . This scale Ir is usually 
larger than the healing length 1^, which plays in this system the role of the Planck length. 
Consequently it is the length Ir which is relevant for the regularization rather than the 



"Planck length" 
Eq. (57) of Ref. 



(healing length). This regularization length can also be deduced from 
However the author did not make it explicit. The regularization 
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removes the divergency in the norm (H8|l and makes the procedure meaningful. The 
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influence of tlie QP can be felt also outside this scale at |x| > Ir since the exponent 7 in 
Eq.f l39|) . contains now corrections due to the QP. They result in a deviation of the high 
frequency tail of the spectrum of Hawking radiation from the Planck black body radiation 
distribution. 

Finally, we want to mention that neglecting QP results in a scale invariance of the 
equation of motion for \k\ — )■ 00. This is the origin of the occurrence of the singular 
solution Xs oc x"'~^ at |a;| — > 0. The regularization we have found close to the horizon 
after accounting for QP is based on breaking this scale invariance for \k\ — t- 00 (see 
Appendix) 
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Appendix A: Function x =^t \x\ — )• 

We try to elucidate here the mathematical features leading to the singularity at the 
horizon and its regularization. Without attempting mathematical rigor we present a 
qualitative discussion. Let us start with the case where the QP is neglected. For a 
stationary solution with variation in the x direction only, we get from Eq. (!32|) 

z/2 - 2iyk 

Ok In Xk 



iak{2i' — 3k — ia) 
for Z/i = or for \k\ 00 

dk\nXk^-lo^ (Al) 
If x{x) is a generalized homogeneous function, i.e. 

x{x) = X^xix), x = Xx 

for all real scaling parameters A, then its Fourier transform fulfills the condition 

Xfc = A(^+^)xfc, k = X-'k. (A2) 
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Making use of (lA2p and = Xdk it is easy to prove that 

~ -7oi 
k 

i.e Eq. (lAip is scale invariant. Taking the QP into account Eq. fl5^ resuhs in 

^InXfc ~ ^/r(U)' (A3) 



from which it follows that 

i.e. Eq. (]A3|) is not scale invariant. This is consistent with the fact that its r.h.s., which 
involves the dimensionless quantity {UkY ■, explicitly depends on the length scale Ir- The 
latter has been found to be the scale below which regularization occurs. Taking in the 
expansion of v{x) and s{x) (close to the horizon) higher order terms in x into account 
higher order derivatives with respect to k will occur in Eqs. ( I30l) . This will prevent the 
decoupling of the fields and Xk and therefore will hamper the analytical calculation of 
^{x,t) and xi^jt)- Nevertheless the breaking of scale invariance still holds and is generic. 
Consequently, the regularization close to the horizon should be robust. 
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